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Abstract 

£T} . Early Mars may have had a warmer and denser atmosphere allowing for the presence of liquid water on the surface. However, 
7— I climate model studies have not been able to reproduce these conditions even with a CO2 atmosphere of several bars. Recent 3D 
< ^> simulations of the early Mars climate show that mean surface temperatures only slightly below 273 K could be reached locally. 

We want to investigate the effect of increased partial pressures of N2 on early Mars' surface temperature by including pressure 

^ broadening of absorption lines and collision-induced N2-N2 absorption. 

AID radiative-convective cloud-free atmospheric model was used to calculate temperature profiles and surface conditions. We 
^^ performed a parameter study varying the N2 partial pressures from to 0.5 bar at CO2 partial pressures between 0.02 bar and 3 bar. 
CNI These values are consistent with existing estimates of the initial, pre-Noachian reservoir for both species. Solar insolation was set 
Cn to be consistent with the late Noachian, i.e. around 3.8 billion years ago. 

^_^ Our ID global mean simulations clearly show that enhanced N2 content in the Martian atmosphere could have increased surface 
Q^ temperatures. An additional greenhouse warming of up to 13 K was found at a high N2 partial pressure of 0.5 bar. Still, even 
PJ at this N2 partial pressure, global mean surface temperatures remained below 273 K, i.e. the freezing point of water. However, 

-h given the magnitude of the N2-induced surface warming and the results of recent 3D studies which show that local mean surface 

O, temperatures are not much lower than 273 K, our results imply that the presence of atmospheric N2 could have led to almost 
1 continously habitable mean surface conditions in some regions. In addition, atmospheric water column amounts increased by up 

^ to a factor of 6 in response to the surface warming, indicating that precipitation might also increase upon increasing N2 partial 
"£j? pressure. 
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1. Introduction 

There is compelling evidence that the early Martian atmo- 
sphere may have been much denser than toda y. For example, 



with pure CO2-H2O atmospheres (e.g., lKastindll99lh . The 
potential warming effect of C O2 clouds has been investigated 



by a number of authors (e.g., [F orget and Pierrehumbertl ll 997 
Mischna et alj|2000l IColaprete and Toonl 120031) . The conclu 



isotope ratio measurements and modelling (e.g.,|McElroy and Yung! sion of these studies was that in order to raise the global mean 
Il976lljakpskv et al.ll994J) and infe rred impact rates from crater 



X 



surface temperature above 273 K, a fractional cloud cover of 

Counting dBrain and Jakoskyl 1998|) suggest intense atmospheric a1mo „ t 100 % with an associated optical depth of the order of 

loss d uring the late Noachian era (see, e.g., review by|jakosky and PhjHips Juld be needed. These conditions are rather unrealistic, 
l2Q0jJ). Outgassing through extensive volcanism might have suggesting that the surface of early Mars would remain at tem- 
peratures below 2 60 K. A recent 3D modelling study of early 
Mars' climate by I Wordsworth et al.l (120131) suggests that even 
for global mean surface temperatures of the order of 250 K 
or lower, locally transient liquid water, alb eit in rather small 



produced a relatively dense CO 2 atmosphere of the order of 0. 1 - 
1 bar (e.g., IPhillips et al.l 1200 ll) . Additional hints for a denser 
early atmosphere come from the existence of jarosite on the sur- 
face (IClark et al. 2005 ) or from extrapolated magma ejection 
trajectories (Manga et al.l 120121) . Together with mineralogical 
(e.g., IPoulet et al. 2005 J) and topog raphical (e.g.. IPerronet al 



amoun ts, may be possible. The same study bv lWordsworth et al. 



2007J or lAnsan and M angold 2006) evidence for liquid water 



on the surface, these observations indicate a probably warmer, 
wetter climate on early Mars before 3.8 Gyr ago. 

Climate model studies of early Mars have encountered im- 
portant problems in reproducing these warm and wet conditions 



d2013l) also showed that, for a 1 bar CO2 atmosphere, annual 
mean surface temperatures of the order of ~260 K are possible 
in some regions. 

An alternat ive to simple CO2-H2O atmospheres was pre 
se nted by, e.g., 
or 



Halevv et al. 



Postawko and Kuhnl fl986l) . I Yung et aD dl997h 
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(120071) . They suggested that SO2 concentra- 
tions of the order of 1 0~ 5 would be sufficient to kee p the surface 
tempe rature above 273 K. Simulations presented bv ljohnson et al 
d2008l) . using a 3D climate model, confirmed these conclusions 
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and implied that S O? could act as an efficient greenhouse gas. 
However, recently, iTian et al.l (120101) suggested that such high 
amounts of SO2 in the atmosphere would lead to aerosol forma- 
tion and subsequent surface cooling. 

In the present Martian atmosphere, N2 is the second most 
abundant constituent. The current N2 partial pressure is about 
2-10 4 bar. Estimates of the total N ? inventory of Mars rang e 
from a few mbars to roughly 0.3 bar (IMcKav and Stokerll989l) . 
Assuming a simple scaling to the atmospheric N2 on Venus re- 
sults in a total N2 inventory of 0.5-0.6 bar for Mars. Two possi- 
bilities have been discussed to explain the fate of the initial in- 
ventory of probably several tens or hundreds of mbars. The first 
involves atmospheric escape. Thermal and non-thermal pro- 
cesses (such as ion pick-up, sputtering etc.) could have eroded 
around 90 % of the initial inventory. By combining all such 
processes, Ijakoskv and Phillipsl (120011) suggest that more than 
99 % of volatiles were lost. Statistical events, such as large 
impacts, could also have significantly eroded the atmosphere. 
The second possibility to explain the lack of atmospheric N2 on 
present Mars, apart from atmospheric escape, is incorporation 
into the crust in the form of, e.g., nitrates. It has been argued 
that several tens of mbar of N2 could hav e been fixed as nitrates 
by, e.g., UV radiation or lightning (e.g., Mancinelli and Baninl 
2003ilSegura and Navarro-Gonzalezl2005l) and thus be removed 
from the atmosphere to the (sub-)surface. The detection of 
such nitrate deposits, if they exist, from orbit is however rather 
difficult, as has been shown by a M ars-analogue study in the 
Atacama desert (ISutter et al.l 120071) . Since this study showed 
that nitrate concentrations increase steeply with depth, in-situ 
drilling missions are probably needed to assess Martian nitrate 
deposits. Combining nitrate data with escape signatures from 
isotopes could then be used to estimate more precisely the ac- 
tual amount of atmospheric N2 in the Noachian period. 

To date, Martian N2 has been mostly of interest in view of 
a potential biosphere and the pos sibility of fixing nitrogen in 



biologically relevant forms (e.g., iBoxe et al.l 120121) . By con- 
trast, the atmospheric implications of high N2 partial pressures 
on the climate of early Mars have not be en investigated so far. 
It has been sh own for (early) Earth (e.g., Ooldblatt et al. 2009, 
Li et al.120091) and the exoplanet candidate GL 58 1 d dvon Paris et 
20101) that increasing N2 partial pressures at fixed CO2 partial 
pressure greatly increases the greenhouse effect, hence surface 
temperatures. Therefore, in this work, we investigate the in- 
fluence of increasing N2 partial pressures on the surface tem- 
perature of early Mars. We want to assess to what extent an 
enhanced N2 content could contribute to warmer surface tem- 
peratures. 

Section [2] briefly describes the model and the simulations 
performed. In Sect. [3] the results are presented and discussed. 
Section|4]presents the conclusions. 

2. Methods 

2.1. Atmospheric model 

We used a one-dimensional radiative- convective model as 
described in Ivon Paris et al.l (120081 I201CH) . The model calcu- 
lates globally, diurnally averaged atmospheric temperature and 



water profiles for cloud-free conditions. In its origin al form, 
the model was introduced by iKasting et al.l dl984allbT) . Model 



updates and improvements to radiative tra nsfer and the tr eat- 
ment of convection are described in, e.g., iKastind dl99ll) . or 
Mischna et alJ (120001) . More details on model equat ions and as- 



sumptions are given in Ivon Paris et al.l (120081 1201 0D 



Model atmospheres are assumed to be composed of CO2, 
N2 and H2O only, and divided into 52 model layers. These 
layers are spaced approximately equidistantly in log (pressure), 
ranging from the surface to the fixed model lid at 6.6xl0~ 5 bar. 
The model calculates the radiative temperature profile by solv- 
ing the radiative transfer equation and imposing radiative equi- 
librium. The radiative fluxes are calculated on two separate 
spectral ranges, one for the stellar radiation (heating rates, mostly 
in the UV, visible and near-IR) and one for the atmospheric and 
planetary surface emission (cooling rates, mostly thermal IR). 
If the calculated radiative temperature gradient is steeper than 
the adiabatic lapse rate, the temperature profile is adjusted to 
the wet adiabatic lapse rate. This wet adiabatic lapse rate in- 
corporates either CO 2 or H2 O as condensing spec ies (see, e.g., 
Kastingll988[lKastin J199ll or lKasting et all 19931 for a detailed 
description of convection in the model). 

The water profile in the model is calculated based on a rel- 
ative humidity (RH) distribution and the water vapor saturation 
pressure. The RH profile is an input to the model (see Sect. 
2.2l for details), whereas the water vapor saturation pressure is 
calculated according to local temperature. Above the cold trap, 
the water profile is set to an isoprofile taken from the cold trap 
value. 

The pressure broadening of CO2 lines is accounted fo r semi- 



empirically by introducing a CO2 foreig n continuum (see | von Paris et al 
' 2OO8U2OIOI) . We follow the approach of lClough et all d 19891) in 
the approximation of line shape and absorption coefficient £ cont 
(in cm 2 ): 



hvc 



Pf 



^contO 7 , T) = v tanh(— -— ) • C fol - e ign(v) — 
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where v is frequency, T temperature, h, c, kg Planck's con- 
stant, speed of light and Boltzmann's constant, respectively, 
alXTfbreign th e continuum (values taken fr om lSchreier and Schimpf 



200 ll and ISchreier and Bottgerll2003l) . and p/, po are the den 



sity of the broadening gas and a reference density, respectively. 
The use of Eq. Q] is consistent with high-resolutio n radia tive 
transfer codes as d eveloped by, e.g., IClough et al.l (1198 lb or 
Clough et al.ldl986l) . For the frequency- and temperature-dependant 



terms in Eq. [T] k distributions on a fixed temperature grid are 
calculated (see 



von Paris et ai1l2008l) 



We also introduced N2-N2 collision-induced absorption (CIA) 
in the IR radiative transfer part of the model. The N2-N2 CIA 
shows two distinct bands, the fundamental band near 4.3 fim 
and the rotational band in the far infrared (A>20 pm). For 
nitrogen-rich atmospheres, this could become a very impor- 
tant contribution to the ove rall IR opacity, hence the green- 
house effect (e.g. on Titan, lLavvas et al.l l2008l IMcKav et alJ 
198 9J). For the fundamen tal band, fits to measurements reported 



by iLaffertv et al.l dl996) were used, whereas for the rotational 
band, absorption coefficients have been taken from molecular 



simulations by Borvsow and Frommholdl (119861) . These stan 



dard ap proximations have been used for example bv lBuehler et al. 



(J2006J) to assess the accurac y of radiative transfer codes com- 
pared to measurements or bv lRinsland et al.l (120101) for retriev- 
ing telluric CO2 profiles based on satellite data. 

2.2. Model scenarios 

All scenarios assumed a solar luminosity S of 0.75 times 
the present value, representing the approximate value 3.8 Gyr 
ago dGoughlll98ll) . 

The model surface albedo was set to As =0.21, which is 



close to the surface albedo of present Mars (e.g., iKieffer et al 



1977) and allows the model to calculate surface temperatures of 
the order of 2 1 7 K when assuming present-day conditions (solar 
insolation, surface pressure). 

The RH profile in the troposphere was set to unity in the 
absence of information on its vertical distribution. This ap- 
proach of assuming a c onstant relative humidity is consist ent 



* 



with other studies (e.g., 'M ischna 

suming RH= 1 and Cola prete and Toonl2003l assuming RH= 



i et all 

Toon: 



2000UKasfinglll991 



as- 

0.5). 



An RH value of unity provides an upper limit for the surface 
temperature since the water vapor content of the atmosphere is 
most likely being over-estimated. 

We varied the partial pressures of CO2 from 0.02-3 bar and 
N2 from 0-0.5 bar. The values chosen for the partial pressures 
are consistent with reported estimates of t he initial Mars volatile 
inventory (e.g. jMcKav and Stokerl 19891) and modeled volca nic 
outgassing throughout the Noachian (e.g. jGrott et al.luOl lb . 



3. Results and Discussion 

3.1. Surface temperatures 

Fig. Q]shows the calculated surface temperatures as a func- 
tion of pco, for the pi\h=0 case. For /?co 2 <2bar, surface tem- 
peratures increase with increasing pco 2 - This is related to the 
enhanced greenhouse effect. However, for pco 2 >2 bar, surface 
temperatures begin to decrease with increasing pco^- This is 
due to enhanced Rayleigh backscattering of solar radiation, which 
favors surface cooling (the planetary albedo increases from about 
0.2 to 0.45, when pco^ increases from 0.02 bar to 3 bar). Hence, 
the increase in greenhouse warming is compensated by the cool- 
ing via increase in Rayleigh scattering, and surface tempera- 
tures decrease upon further increasing pco 2 - 

Also plotted in Fig. [TJare annual mean surf ace temperatures 
from 3 D simulations of early Mars presented by IWords worth et al. 
(120131) . Surface temperatures agree rather well over the consid- 
ered range of CO2 partial pressures. Note, ho wever, that we 
use an a pproxima tion of the CO2-CO2 CIA from lKasting et al. 
dl984bl) . whereas IWordsworth et al. ( 20 1 3 ) use a more recent 
parameterization from lWordsworth et al.l(l2010l) . It was demon- 
strated by that study that this new parameterization might lead 
to a cooler surface due to less CO2 opacity. Furthermore, our 
model is a cloud-free code, whereas IWordsworth et al.l (120131) 
include CO2 clouds. In their study, these lead to an increase 
in surface temperature. Fig. [TJ suggests that these two effects 
approximately cancel. 
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Figure 1 : Comparison of calcu lated surface temp e ratures (at pn 2 =0) with an- 
nual mean temperatures from IWordsworth et al.l J2013I) (their Fig. 2, blue 
points). 
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Figure 2: Surface warming in K with respect to pn 2 =0 bar scenarios. 

Fig. [2] shows the surface warming as a function of p^ and 
/?cOo with respect to the respective p^-Q scenario. It is clearly 
seen that all scenarios are warmer than the p^-Q case. 

The increase of p^ has a significant effect on surface tem- 
peratures for CO2 partial pressures of 0. 1 bar and higher, be- 
ing most pronounced at CO2 pressures of the order of hun- 
dreds of mbar. The maximum effect upon increasing p^ 2 to 
0.5 bar is nearly 10 K (at pco^ -2 bar) and 8.5 K(atpcOo=3bar). 
When increasing p^ 2 to 0.3 bar (which constit utes the upper 
limit from estimates in lMcKav and Stokerl 19891) . the additional 
greenhouse effect already results in 6 and 7 K higher surface 
temperatures, respectively. At 0.5 bar of CO2, the effect on sur- 
face temperature of an increase in N2 is even more substantial, 
namely 9 K at p^-, =0.3 bar and 13 K at /?n 2 =0.5 bar. 

An increase of p^ 2 leads to an increase of Rayleigh scat- 
tering of incoming solar radiation, hence an increase in plane- 
tary albedo and an associated negative radiative forcing. How- 



ever, since CO? is a much m ore efficient scatterer than N2 (e.g., 
Vardavas and Carverll984l) . the radiative forcing due to Rayleigh 
scattering is rather weak at high pco, ■ The maximum radiative 
forcing in the parameter range considered occurs at /?n, =0.5 bar 
and decreases from -7.4 Wm -2 at pco, =0.02 bar to -1.1 WrrT 2 
at pco^=3bar. Due to higher surface temperatures upon in- 
creasing pn 2 , the H2O content of the atmosphere increases (see 
below), which leads to increased near-IR absorption of solar 
radiation in the H2O bands. This effect introduces a positive ra- 
diative forcing and partly compensates (by about 10-20%) the 
negative forcing due to Rayleigh scattering. 

At the surface, the increasing p^ leads to a positive ra- 
diative forcing for thermal radiation, i.e. an increasing green- 
house effect. This radiative forcing is maximal at /? N2 =0.5bar 
and also decreases with increasing /?co,, as does the Rayleigh 
forcing. At pcoi=0.5bar, it amounts to 11.6 WrrT 2 , whereas 
at pco 2 =3bar it is only 2Wm 4 . This behaviour is due to the 
fact that at CO2 partial pressures above roughly 1 bar, the near- 
surface atmospheric layers are already optically thick at all IR 
wavelengths, even without the presence of atmospheric N2. There- 
fore, the surface warming shows a maximum at pco 2 =0.5-1 bar. 
The N2-N2 CIA (see Sect. [27Tb contributed less than 1 K to 
the increase in surface tem perature. Pressure broadening of ab- 
sorption lines by N 2 (e.g., boldblatt et alj |2009l 10 et al.ll2009l 
von Paris et al.ll2.Q10l) constitutes the far more important effect. 
Further sensitivity tests which varied the amount of H2O and the 
strength of the CO2-CO2 CIA indicated that the warming effect 
of N2 is nearly independent of other model radiative transfer 
details. 

Therefore, our calculations suggest that N2 could play an 
important role for warming early Mars by providing a signifi- 
cant additional greenhouse effect. 



3.2. Implications for liquid surface water and habitability 

Geological and morphological evidence regarding valley net- 
works on Mars implies that during the Noa chian, early Mars 
had a period of extensive fluvial activity (e.g., Pou et et al.l2005 



Ansan and Mangoldl2006llFassett and HeadOOl II) . Climate mod- 
eling studies have shown that liquid water could exist temporar- 
ily and locally on early Mars even at very low ann ual mean sur- 
face temperatures (e.g., IWordsworth et al.ll2013l ). The annual 
amounts of liquid surface water, however, are of the order of 
0.1-1 kg m -2 , i.e. far less than needed to explain the observed 
fluvial features. 

Fig. [3]shows the atmospheric H2O column, as a function of 
/?C0t f or different p^ . For reference, we show also the value of 
the mean atmospheric H2O column on modern Earth as a hori- 
zontal line. Note that on modern Earth (1990), /?co 2 =3-5T0~ 4 bar. 
It is clearly seen that an increase of N2 partial pressure leads 
to considerably wetter atmospheres, as could be expected from 
Fig. [2] The increase of atmospheric H2O reaches up to a factor 
of 4-5 (at 1, 2 and 3 bar of CO2) when increasing N2 to 0.5 bar, 
and around a factor of 6 at 0.2 bar of CO2. Therefore, it seems 
likely that precipitation and annual surface liquid water could 
also increase by corresponding amounts. Note, however, that 
even the considerably increased H2O columns at p^ =0.5 bar 
only represent about 10~ 3 -10~' of the terrestrial atmospheric 




Figure 3: Atmospheric H2O column as a function of CO9 and N2 partial pres- 
sures. Value of Earth mean column indicated by horizontal line. 



H2O column. Nevertheless, this increase could help to resolve 
partly the apparent contradiction between observed early fluvial 
activity and climate modeling studies. 

Also, the calculated increase in global mean surface temper- 
ature of up to 13 K for significant amounts of atmospheric N2 
implies that locally, mean annual temperatures coul d be close to 
or hig her than 273 K (see, e.g., temperature maps in lWordsworth et al 
2013b . These surface temperatures would allow for liquid wa- 
ter, one of the key requisites for life on Earth. The presence of 
liquid water is also a fundamental part of virtually all habitabil- 
ity definitions (e.g., Kasting et alJ 1993 fo r terre strial planets, 
Mars Exploration Program Analysis Groupll2005l for the Mars 
exploration context). Therefore, our results imply that early 
Mars could have been continously habitable at least in some 
regions. Due to the (chaot ic) evolution of bot h Martian eccen- 



tricity and obliquity (e.g., lLaskar et al.l 120041) . the location of 



these continously habitable regions might change on timescales 
of 10 4 -10 5 years (see, e.g.. IWordsworth et al.ll2013l) . Further- 



more, it should be noted that with increasing surface pressure, 
the dependence of surface te mperature on topography would 
most likely also increase (e.g.. IWordsworth et al.ll2013l) . hence 
complicating the explanation of the mostly high-altitude late 
Noachian valley networks. 

3.3. Noachian atmospheric N2 content 

Our parameter study aimed at investigating the influence of 
atmospheric N2 on surface temperatures. Fig. |2] shows that 
significant amounts of N2 of the order of 0.3 bar are needed in 
the model to warm effectively the surface. Such amounts are to- 
wards the upper limit of Martian N2 estimates (iMcKav and Stokei 
1989). Simulations of atmospheric escape suggest that most of 
the early protoatmosphere was probably lost very shortly af- 



ter th e formation of Mars (e.g., iTian et alj|2009l Lammer et al 



20131 and references therein). The subsequent re-formation of 

a secondary atmosphere containing significant amou nts of N2 

could be possible by, e.g., impact-related delivery (e. g.. Manning et al 



20081 l2009[lde Niem et al.l2012h . since chondrit e-type impactor s 



contain relatively high amount s of nitrogen (e.g.. lKerridgell985 
Hashizume and Sugiuralll995l) . Volcanic outgassing or nitrate 



decomposition has also been sugge sted to contribute to sec 
ondary N2 in the atmosphere (e.g. , lOwen et al.l 1 19771 IWallis 



19891 IPemnll 1 994 Ijakoskv et al.ll 1994 

Hence, upcoming investigations of the Martian surface and 
sub-surface (either from orbit or in-situ) should ideally include 
the search for potentially missing N2. Future atmospheric es- 
cape or outgassing studies should also consider N2. In doing 
so, the present-day Martian N2 reservoir could be constrained 
much better and it could be assessed whether N2 was present 
in high enough concentrations in the early Mars atmosphere to 
influence the Noachian climate. 

4. Conclusions 

We have presented model calculations of the early Mars at- 
mosphere with a ID cloud-free atmospheric model. We varied 
N2 partial pressures as well as CO2. To our knowledge, ours 
is the first study in which the radiative effects of N2 (pressure 
broadening and N2-N2 collision-induced absorption) were in- 
vestigated in simulations of the early Martian climate. 

Results suggest that increasing N2 partial pressures could 
lead to significant surface warming, without accounting for CO2 
clouds or additional greenhouse gases. Surface temperatures 
increased by up to 13 K at high N2 partial pressures of 0.5 bar. 
Nevertheless, global mean surface temperatures still remain be- 
low 273 K. Given the results of previous 3D studies of CO2- 
H2O atmospheres, however, our results suggest that local an- 
nual mean surface temperatures above the freezing point of wa- 
ter might be possible. The water column density increased by 
up to a factor of 6, indicating that precipitation might also sig- 
nificantly increase. 

Detailed studies of the outgassing history, impact delivery 
and atmospheric escape of N2 are needed to assess whether high 
N2 partial pressures could be consistent with the late Noachian 
period. 

A further investigation of the effect of increasing N2 upon 
the transient existence of liquid water on the surface and the for- 
mation of fluvial features with detailed 3D atmospheric models 
is warranted in the future. 
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